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r^,' Abstract 

We analyze neutrino oscillation for the general case when the initial neu- 
■ trino is not in a pure flavor state. We show that, after such a neutrino beam 

, propagates for a while, the probability of detecting any pure flavor state de- 

pends even on the CP-violating Majorana phases in the mixing matrix. The 
I/"") ■ dependence remains even when energy spectrum of the initial beam is taken 

. into account. We discuss various implications of this dependence. 



Through solar and atmospheric neutrino data as well as ground-based oscillation 
experiments, we now know that neutrinos have mass, and they mix [1]. Two of the 
three mixing angles are already known to be non-zero, and in fact quite large, 
whereas for the other there exists only an upper bound. The masses themselves are 
not known, and cannot be known from oscillation data alone, although we have a 
good idea of the mass differences (more precisely, the differences of mass-squared 
values) between different eigenstates. 

Mass matrices contain information not only about mass eigenvalues and mixing 
angles, but also about phases responsible for CP violation. The number of such 
phases depend on whether the neutrinos are Dirac fermions or Majorana fermions. 
For N generations of leptons, there are ~(N — 1)(N — 2) Dirac phases. In case of 
Majorana neutrinos, the number of phases is ^N(N-l), i.e., N—l more phases com- 
pared to the Dirac case [2-4] . These extra CP- violating phases are sometimes called 
"Majorana phases" , and we will use this terminology for the sake of convenience. 

An intriguing question is how to observe these phases should they exist, i.e., 
if the neutrinos are Majorana particles. Since Majorana neutrinos are obtained in 
theories where lepton number is violated, the observability of the Majorana phases 
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was discussed in the context of lepton number violating processes [4, 5] , and for a 
long time it was believed that these phases can be observed only in such processes. 
Later, it was shown [6,7] that lepton number violation in the process is not necessary 
in obtaining information about the Majorana phases. 

Inspired by this knowledge, we can ask the question whether the Majorana phases 
can be observed in neutrino oscillation experiments. This is what we try to do in 
this article, for neutrinos oscillating in the vacuum and also in matter. 

It is best to address this question assuming there are just two generations of 
neutrinos, where the analysis can be carried out analytically. In this case, the 
mixing matrix contains only one phase and it is of the Majorana type. The flavor 
states v e and are related to the mass eigenstates through the relation 

/ v e \ _ / cos 9 e ia sin 9 \ / V\ 

\ v ii) v — e ~ ia s ^ n ^ cos ^ / V V2 

Suppose now we start with a monochromatic neutrino beam, with energy E, which 
is not necessarily in one of the flavored states, but is a general superposition of the 
two flavors: 

|^(0)) = A\u e ) + B\u,) (2) 
= {Acos9 -Be- ia sm9) \v x ) + (Ae ia sin9 + B cos9) \u 2 ) , (3) 

where A and B are real, and A 2 + B 2 = 1. We could have been more general and 
assumed that there is a relative phase between A and B, but that is not necessary 
for the argument that we are going to present. Time evolution of this state gives, 
apart from an overall phase that is unimportant, the result 

\ip(t)) = (Acos9 - Be- ia sin9) \v x ) + (Ae ia sin9 + B cos9)e- 2iS \v 2 ) , (4) 

where 

2 2 

mi — mi 
<5 = „^ t, 
AE 

m\ and m 2 being the mass eigenvalues of v\ and v 2 

We can now ask the question what is the probability of finding one of the flavors, 
say u e , in this resulting beam. The answer is given by 

P Ve {t) = \(Ve\m)\ 2 ■ (6) 

Using Eq. (pQ) to express v e in terms of the eigenstates and performing a few steps 
of trivial algebra, we obtain 

P Ue {t) = A 2 {1- sin 2 29 sin 2 5) + B 2 sin 2 29 sin 2 5 

-2AB sin 29 sin 5 (cos 2 9 sin(a + 5) + sin 2 9 sin(a - 5)) . (7) 

Clearly, if A = 1, this will denote the survival probability of z/ e 's in a beam that 
started as z/ e , and in this case our expression reduces to the familiar expression for 
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that case. Our result also agrees with the analysis of neutrino oscillation probabilities 
performed [8] for A and B both non-zero but assuming Dirac neutrinos, which 
implied a = 0. The same analysis can hold for Majorana neutrinos if CP violation 
is neglected. But our expression is more general than all these results, and it shows 
that if neither A nor B vanishes, i.e., if the initial beam is not in a pure flavor state, 
then the oscillation probability depends on the Majorana phase a. 

What is more important, the dependence on a does not get washed out even if 
we have a neutrino beam with a large spread of energy which travels through a long 
enough duration of time. In this case, we need to take averages of the 5-dependent 
quantities over the energy spectrum. Quantities like sin 5 or cos 5 vanish under such 
averaging, but sin 2 S averages to |, and we obtain 

P Ue {t) = A 2 {I - - sin 2 26) + ^B 2 sin 2 20 - AB sin 26 cos 26 cos a . (8) 

We now proceed to derive the corresponding formulas for oscillations in matter. 
We denote the 2x2 Hamiltonian in the flavor basis by 

\tt 2 i ri 2 2 J 

The Hamiltonian has to be hermitian, irrespective of whether we have Dirac type of 
neutrinos or Majorana type, and so we must have Hn and H22 real and H21 = Hl 2 . 
Now, if there is a matrix U of the form shown in Eq. ([1]) that makes U^HU diagonal, 
the parameters of this matrix are given, in terms of the elements of H, by 

6 = l - tan" 1 ( ff^f^J > « = arg(#i 2 ) • (10) 

For neutrinos passing through matter, the Hamiltonian depends on the density [9]. 
Density corrections appear in the diagonal elements of the Hamiltonian written in 
the flavor basis. As a result, a does not change with density, although 6 does. 

So now, suppose the initial state produced is the superposition of v e and u^, as 
given in Eq. (j2J). At the production point, the density is such that the effective 
mixing angle is 6 . Thus, 



cos 6*o e ta sin 9 \ / V\ 
Up J \ —e~ ia sin 6 cos 6 J \ v 2 



where the subscript "0" indicates the production point in matter. Thus, 

|V(0)) = (Acos6 Q - Be- ia sin6 ) \u 1 } + (Ae ia sin6 + B cos6 ) \u 2 } . (12) 

The probability that this state is in the first eigenstate is given by |v4cos^o — 
B sin 6oe~ ta \ 2 . In the adiabatic case, this probability remains constant, although 
the eigenstate, as a superposition of the flavor states, changes with density. If at 
the end of the journey the neutrino is detected in the vacuum where Eq. (CQ) holds, 
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the probability of finding a \v e ) in this state is just cos 2 8. Similarly adding the 
contribution from the other eigenstate, we obtain 

P { jf\t) = \Acos6 - Bsm6 e- ia \ 2 cos 2 8+ \Asin8 e ia + B cosd \ 2 sin 2 8 

= -A 2 (1 + cos 28 cos 28) + -B 2 (1 - cos 28 cos 28) 

-AB sin 20 o cos 26 cos a . (13) 

The superscript 'ad' reminds us that this is the result if the adiabatic condition, 
described above, holds. This is, of course, the energy-averaged expression, because 
we have neglected the oscillatory terms by working with an explanation in terms 
of probabilities [10] and not of amplitudes. But even in this expression, there is an 
a-dependence. 

If the conditions are non-adiabatic, then neutrinos can jump from one eigenstate 
to another. If this jump probability is called Pj U mp, then the probability of finding 
u e would be given by [11, 12] 

Pu. = (i - iW4 ad) + Wi - Pt ] ) ■ (14) 

The jumping probability depends on the expression of the eigenvalues in terms of the 
parameters that appear in the Hamiltonian in the flavor basis. These eigenvalues are 
independent of a and therefore the jumping probability is independent of a as well. 
However, the expression for P Ve still depends on a through the adiabatic probability 
that appears in this expression. 

There are claims in the literature that the Majorana phases are unobservable in 
neutrino oscillation experiments, whether in vacuum [2-4] or in matter [13]. How- 
ever, these claims have all been made in the context of initial neutrino states of pure 
flavor. This is consistent with our results since the a-dependence in all relevant 
formulas appear with a co-efficient AB, and therefore vanish for pure flavor states 
which have either A = or B = 0. 

Presence of a non-zero value of the CP-violating phase a can introduce some 
qualitative features in neutrino oscillation probabilities. As an example, suppose we 
have an initial neutrino beam with some known admixture of two flavors. If a = 0, 
the probability of finding v e after this beam travels for a time t in the vacuum is 
given by 

P (t) = A 2 - [( A 2 - B 2 ) sin 2 28 + 2AB sin 28 cos 2d] sin 2 5. (15) 

If the initial state was such that both A and B are positive with A > B, the right 
hand side of this expression is always less than A 2 , assuming 8 < 71/ A which is 
certainly valid for v e -v^ oscillation, as inferred from various solar and terrestrial 
experiments [1]. However, if a 7^ 0, we need to use the expression given in Eq. 
(EJ), and P Ve {t) is no longer guaranteed to be less than A 2 . Thus, if one observes a 
reinforcement of the dominant component of the beam in an oscillation experiment, 
it definitely signals a non-zero value of the phase a. 

Let us conclude this article with a qualitative feeling for the magnitude of the 
importance of a oscillation formulas, which might help in devising experimental 
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Figure 1: Various equal probability contours. Details have been given in the text. 



techniques for finding a from neutrino oscillation experiments. We will do this by 
considering fictitious data coming from two different oscillation experiments. It has 
been argued [14] that the analysis of data from a single oscillation experiment is 
conveniently done by introducing the dimensionless parameter 



2 2 
777.2 — m l 



^^eT*' (16) 

where (E) is the average energy in the incoming beam. Different experiments can 
have different average energy in the incoming beam and different distances between 
the source and the detector, and consequently different values of 5 even when they 
are measuring oscillation between the same two flavors. In order to provide an 
illustrative example of the point we are discussing here, we consider, for the sake of 
simplicity, that two experiments have the same (E) and the same source-detector 
distance (or at least the same ratio of the two quantities just mentioned), and 
therefore the same value of 5. One of these experiments is working with an initial 
beam that is purely z/ e , and observes that P Ue = 0.8 at the detection point. We can 
ask which values of 5 and the mixing angle can make it possible. The answer will 
be independent of a since the initial beam was purely v e , and has been shown in 
Fig. [TJ The shape of this line of course depends on the energy spectrum of the initial 
beam. For the sake of definiteness, we have taken a Gaussian energy spectrum with 
a standard deviation equal to 0.2(E). 

Now suppose that in the second experiment, the initial neutrino beam has A = 
0.8 and B = 0.6 in the initial beam, and at the detector, one finds P Ue = 0.6. The 
energy spectrum is the same. If a = 0, the equal probability contour is given by 
the line that reaches leftmost in Fig. [TJ It is worth pointing out a special feature of 
this line compared to the line obtained from the first experiment. Near the bottom 
end of the line, we have two solutions of the mixing angle for the same value of 5. 
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This happens because, when A and B are both non-zero, the ^-dependence of the 
expression in Eq. (j7J) is not monotonic. 

Anyway, this line is clearly inconsistent with the result of the first experiment. 
However, if a 7^ 0, this need not be the case. We show the contour for a = it/ 2 in 
Fig. [T] with a thick solid line, which has clearly some intersections with the curve 
from the first experiment. Thus, non-zero values of a will make the results of the two 
experiments consistent with each other, and the value of a can even be determined 
from such data from two experiments. 

Real experiments, of course, will have some error bar on the detected probability, 
so there will be a band rather than a contour corresponding to a single experiment. 
There will be other complications in the analysis because the energy spectra of no 
two experiments will be the same. Obtaining an initial beam which is not a flavor 
eigenstate is not straight forward as well. Here, we have only discussed the matters 
of principle and indicated possible ways in which the CP-violating phase a might 
be detected from neutrino oscillation data. 
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